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Abstract. It is generally accepted that relativistic jet outflows power the nonthermal emission from active galactic 
nuclei ( AGN) . The composition of these jets - leptonic versus hadronic - is still under debate. We investigate the 
microphysical details of the conversion process of the kinetic energy in collimated relativistic pair outflows into ra- 
diation through interactions with the ambient interstellar medium. Viewed from the coordinate system comoving 
with the pair outflow, the interstellar protons and electrons represent a proton-electron beam propagating with 
relativistic speed in the pair plasma. We demonstrate that the beam excites both electrostatic and low-frequency 
magnetohydrodynamic Alfven-type waves via a two-stream instability in the pair background plasma, and we 
calculate the time evolution of the distribution functions of the beam particles and the generated plasma wave 
turbulence power spectra. For standard AGN jet outflow and environment parameters we show that the initial 
beam distributions of interstellar protons and electrons quickly relax to plateau-distributions in parallel momen- 
tum, transferring thereby one-half of the initial energy density of the beam particles to electric field fluctuations 
of the generated electrostatic turbulence. On considerably longer time scales, the plateaued interstellar electrons 
and protons will isotropise by their self-generated transverse turbulence and thus be picked-up in the outflow 
pair plasma. These longer time scales are also characteristic for the development of transverse hydromagnetic 
turbulence from the plateaued electrons and protons. This hydromagnetic turbulence upstream and downstream 
is crucial for diffusive shock acceleration to operate at external or internal shocks associated with pair outflows. 
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1. Introduction 

The detection of intense medium-energy gamma radia- 
tion from over 60 blazar active gactic nuclei (hereafter 
abbreviated as AGNs) with the EGRET instr ument on 
the Compton observatory (Hartman et al. 1999| ) an d TcV 
gamma radiation from several BL-Lac AGNs (Pohl [2001 ) 
shows that nonthermal gamma-ray production is a signif- 
icant dissipation mechanism of jet energy generated by 
black-hole accretion. Besides the modelling of the broad- 
band nonthermal radiation, the composition of the jet 
plasma - i.e. electron-positron pair jets (leptonic jets) ver- 
sus electron-proton jets (hadronic jets) - and the acceler- 
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ation of these jet particles to relativistic energies are main 
subjects of the current theoretical research. 

Gamma radiation in leptonic models of broadband 
blazar emission is attributed to synchrotron self-Compton 



(Maraschi et al. 1992; Bloom & Marscher 1996; Tavecchio 



et al. 1998) or external Compton (Dermer & Schlickeiser 
|1993| Sikora et al. |1994| ; Bottc her et al. |1997| ; Dermer 
et al. 1997; Arbeiter et al. |2002D p rocesses (see Bottcher 
[2001 and Sikora & Madejski 2001 for recent reviews). In 
hadronic models, secondary photopairs and photomesons 
and secondary mesons are produced when energetic pro- 
tons and ions interact either wi th am bient synch rotro n 
photons (Mannheim & Biermann 1992[ Mannheim 1993 ), 
photons of the exter nal fie ld (Bednarek & Protheroe 1999 ; 
Atoyan & Dermer [2000 ) and/or ambient matter fields 
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(Pohl & Schlickeiser 2000| ). Observationally, future detec- 
tion of high-energy neutrino emission correlated with high- 
energy photon emission (Schuster et al. 2002) will provide 
the ultimate test between the lcptonic and hadronic jet 
models. 

Most existing radiation models of AGN jets are very 
unspecific on the microphysical details of the conversion 
of the kinetic jet energy into energetic charged particles 
and subsequently into radiation. Without detailed discus- 
sion these models often assume that a significant fraction 
of the accreted kinetic energy is injected into nonthermal 
pairs and/or hadrons with power-law distribution func- 
tions. This efficient conversion is attributed to the scenario 
that the outflowing relativistic jet plasma has produced a 
relativistic shock with fully developed hydromagnetic tur- 
bulence in order to allow for efficient non-thermal diffusive 
particle acceleration at the collisionless shock fronts: 

(1) This neglects the fact that it takes a finite time to 
build up the necessary turbulence in the two-stream multi- 
fluid system consisting of the relativistically moving jet 
plasma and the traversed interstellar or intergalactic hy- 
drogen plasma. 

(2) ft is not clear from the beginning that the out- 
flowing relativistic jet plasma will generate primarily 
(and enough) transverse magnetohydrodynamic turbu- 
lence, which is crucial for efficient particle deflections in 
the up- and downstream region of the shock waves. It is 
well conceivable that most of the kinetic blast wave energy 
is transferred to electrostatic plasma turbulence and not 
to transverse hydromagnetic turbulence. 

(3) The properties of the generated plasma turbulence are 
decisive both for the formation and the nature of the de- 
veloping collisionless shock waves. It is known from non- 
relativistic shock theory that the inclusion of the finite 
pressure and energy flux of the generated plasma turbu- 
lence in the Rankinc-Hugoniot shock relations strongly 
modifies the standard fluid shock pro pertie s (Vainio & 
Schlickeiser |1999| , |200l| ; Lerche et al. [2000D and subse- 



quently the energy spectrum of the accelerated nonther- 
mal particles. 

It is the purpose of this work to consider more thor- 
oughly some of the microphysical details of the energy 
conversion in relativistic jet outflows. We consider the 
energisation of relativistic particles in the jet by interac- 
tions with the surrounding medium following the earlier 
work of three of us (Pohl & Schlickeiser 2000; Pohl et al. 
2002). There the AGN jet has been assumed to be a cloud 
of dense electron-proton plasma which moves relativisti- 
cally through the electron-proton interstellar medium of 
the AGN host galaxy. The plasma cloud is assumed to 
have a cylindrical shape with thickness d* 
which is small compared to the radius r* = 1 
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the cylinder (see Fig. I of Pohl & Schlickeiser (2000) for a 
sketch of the assumed cloud geometry) . It was shown that 
in such hadronic jets swept-up ambient matter is quickly 
isotropised in the jet cloud frame by a relativistic two- 
stream instability, which provides relativistic particles in 
the jet cloud without invoking any acceleration process. 



Here and in the following the index * indicates the val- 
ues of physical quantities in the laboratory (AGN) frame; 
quantities not indexed are in the jet frame. 

Here we study the leptonic variant of this jet out- 
flow model, i.e. we model the outflowing jet cloud as a 
one-dimensional channeled outflow of thickness d* with 
relativistic bulk velocity V, consisting of pairs of elec- 
trons and positrons of density n\ instead of electrons and 
protons. To avoid dramatic pair annihilation, the den- 
sity of pairs must be limited, and we allow for a ther- 
mal pair distribution with non-relativistic temperature 







fcBT pa i r /(m e c ) << I in the jet rest frame. This 



beam of pairs propagates into the surrounding interstellar 
medium that consists of cold protons and electrons at rest 
of density n*. 

For mathematical convenience we assume that the out- 
flow is directed parallel to a uniform background magnetic 
field. The assumption of the magnetic field directed along 
the ejecta velocity enormously facilitates the analytical 
treatment of the two-stream instability, but crucially de- 
pends on the location of the particle energisation with 
respect to the large-scale magnetic field and the confine- 
ment of the jet ejecta. In magnetohydrodynamic models of 
jets in accreting systems the poloidal magnetic field com- 
ponent is more strongly (oc R~ 2 ) suppressed by the side 
expansion of the ejecta than the transverse magnetic field 
component (oc -R -1 ), so that our assumption will hold at 
close distances to the central object. Moreover, in the pres- 
ence of a dominating transverse magnetic field it can be 
argued that the ambient medium will penetrate the ejecta 
by one Larmor radius only and the momentum exchange 
should be faster and more efficient than in the case of a 
dominating poloidal magnetic field considered here. ^ 

Viewed from the coordinate system comoving with the 
pair outflow, the interstellar protons and electrons repre- 
sent a proton-electron beam propagating with relativistic 
speed — V antiparallel to the uniform magnetic field direc- 



tion. Modifying the analysis of Pohl & Schlickeiser (2000) 
and Pohl et al. ( 2002| ) for cold electron-proton outflows to 
finite temperature pair outflows, we demonstrate that the 
beam excites both electrostatic and low-frequency magne- 
tohydrodynamic Alfven-type waves via a two-stream in- 
stability in the pair background plasma. We study the 
time evolution of the beam particles, the generated plasma 
wave turbulence power spectra and discuss the radiation 
signatures of such systems. 



2. Basic equations 

2.1. Outflow parameters 

In the comoving frame the total phase space distribution 
function of the plasma in the blast wave region at the start 
(time t = 0) is 

f(p,t = 0) = mfi(p,iJ,,t = 0) + n b f b (p,n,t = 0) = 



We are grateful to the referee for noting this difference. 
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mS((i + i)s( P -p) + n b T( P ,e) 



2tt P 2 

with the thermal pair distribution 

T(p,Q) = [4Tr{m e c) 3 eK 2 (Q~ 1 )]~ 1 exp(- 
and the dimensionless pair temperature 



1 + 

7TI-I-- 



e 



(i) 



,(2) 



(3) 



Here /i = pu jp is the cosine of the pitch-angle of the parti- 
cles in the magnetic field B, and P — TmV — mc\fT 2 — 1 
where V = l/yjl — (V 2 /c 2 ). The number densities trans- 
form as rii = Tn* , n b = n* b /T. 

Earlier radiation modelling of AGN jets (e.g. Dermer 
& Schlickeiser [19 9$ Pohl & Schlickeiser [2000fl indicated 
that the following parameters are appropriate: the density 
of pairs in the jet cloud is of order nj* 

18 
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implying n b = W b n b _s = 10 8 n^ iqT^ 1 cm ~ 3 , and much 
larger than the interstellar gas density n*. The initial 
Lorentz factor is of order V = 10 2 r 2 and the total number 
of pairs in the plasmoid is N tot = 3 • 10 51 (r 14 ) 2 dl 3 n* b 10 , 
corresponding to a total kinetic energy of E tot = 5 • 
W i7 T 2 (r* u ) 2 dl 3 n* bw erg. 

The background magnetic field in the blast wave 
is parametrized by the equipartition parameter es = 
0.1es,_i through 



B = 0.14(e B! _in M e_ 4 ) 1/2 Gauss 



(4) 



pair plasma 
g p = 7.3 • 10 



-4 1/2^, 



-3/2 



1.6 • lV-\nt w yi*T- 2 l '*Q-_r « 1 



(5) 



is much smaller than unity, so that the kinetic plasma 
description is appropriate. 

The beam distribution function /j in Eq. (Q) is unsta- 
ble and excites electrostatic and low-frequency transverse 
plasma waves. We want to calculate the time scales t e 
and t t it takes these plasma waves to relax the incoming 
interstellar proton-electron beam to a stable distribution 
with respect to both instabilities. Of particular interest is 
the timescale t t for beam isotropisation: if this relaxation 
time is much smaller than the light travel time d/c then 
an isotropic distribution of protons and electrons in the 
blast wave region is efficiently generated. 

2.2. Parallel plasma waves in hot pair plasmas 

We restrict our analysis here to parallel propagating 
electrostatic and low-frequency transverse plasma waves. 
Because 



m/n b = T 2 n*/n* b = 10 6 T 2 n*/nt 10 << 1 , 



(G) 



the beam is weak. Therefore the contributions from the 
beam to the plasma wave dispersion relations of elec- 
trostatic waves and transverse waves at frequencies ujr 
well below the non-relativistic electron gyrofrequency 
{\ojr\ << |^ e |) arc perturbations to the longitudinal 
(5iA e = 0) and transverse dispersion relations (SftAi = 0) 
in the thermal pair plasma. The properties of the longitu- 
dinal and transverse dispersion relations in a thermal pair 
plasma are derived in Appendix A. 

From the real part of the longitudinal dispersion rela- 
tion in the thermal pair plasma we obtain the electrostatic 
waves 



(7) 



in the wavenumber range \k\ < 2u p , e /(c@ 1 / 2 ) as the only 
longitudinal plasma mode. These waves are not damped 
(luj = 0) in the super luminal wavenumber range |fc| < 
v^ojp.e/c, whereas in the subluminal wavenumber range 
\[2uj p ^jc < \k\ < 2ui Pi e/(cQ 1 ' 2 ) the damping rate is 



Ui = -7rsgn(wR,) 



ck[c 2 k 2 



2c^ c ]eK 2 (i/e) 



x exp 



c|fc| 



Ojc 2 k 2 -2u;l e 



29 
c\k\ 



c 2 k 2 - 2u 2 



29(c 2 fc 2 - 2io 2 p e ) 



2 2 k 2 



(8) 



for a non-relativistic initial normalised (co-moving) pair where 
temperature = 10 0_4. The plasma parameter of the 



w p , e = ^4Tre 2 n b /m e = 5.64 • 10 8 K 10 ) 1/2 r 2 1/2 Hz (9) 

is the electron plasma frequency. 

Likewise, from the real part of the transverse disper- 
sion relation in the thermal pair plasma we obtain the 
transverse dispersion relation 



2 v 2 n 2 e k 2 
ujR ~ ft 2 + v 2 k 2 

where 

fl e = — = 1.76 • 10 7 (B/Gauss) Hz 

m e c 

= 2.46 • 10 6 (e B ,-inb,80-4) 1/2 Hz 



(10) 



(11) 



denotes the absolute value of the nonrelativistic electron 
gyrofrequency and 



V e = 



\/2uj p 



9.2 ■ 10 7 (e B ,_ie_ 4 ) 1/2 cms" 



(12) 



is the Alfven speed in the pair plasma. In Eq. (|10[ ) we 
adopt the convention that positive values of lor > denote 
left-handed circularly polarized waves, and that negative 
values of lur < denote right-handed circularly polarized 
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For small wavenumbers \k\ << k e , where Integrating Eq. (^) over time gives 

k e ee Q e /V e = 2.7 • lCT 2 ^ 2 cm- 1 (13) j e(M ) _ J e ( M = 0) = 2tt 2 V u, 2 , / dp,, / 

is the inverse pair skin length, Eq. (nO) reduces to 

K 2 fc 2 (14) 



P U Pj -§( fc ^ll 



7 



7m, 



o 



dp,. 



(19) 



whereas in the opposite case |fc| >> k e Eq. ( |l0| ) ap- 
proaches the pair-cyclotron wave limit ~ J7 2 . For small 
wavenumbers (|fc| << k e ) we obtain for the pair Alfven 
waves the damping rate 



At the same time the phase space density of the beam 
particles evolves as 



.,. tt 1 / 2 n e v e k e 

UJl,A{k) ~ TTT ex P 

2 V t h,e \k\ 



Vth,ek 



dfi 

at 



dpn 



dk I(k,t)S(to R -^) dfl 



(15) 



where vth,e — v20c = y / 2kBT e /m e is the non-relativistic f, 
thermal pair velocity. 



jrrii dp\\ 

In terms of the phase space distribution function 
1 



27rpj_ 



(20) 



(21) 



2.3. Time evolution of the electrostatic instability 

Neglecting spatial dependencies, the time-dependent be- 
haviour of the intensities of the excited electrostatic waves 
is given by 



the two coupled Eqs. (|T^) and (^(J) read 

I e (k,t)-I e (k,t = 0)=*J2 ulif 
kp\\ 



dpn 



l + (-5L)2 



dl e (k,t) 
dt 



2ip e I e (k,t), 



(16) 



where the growth rate ipe is (Schlickeiser [2002 , Ch. 11.2) 

27T 2 



and 



l + (-£L)2 



1 dt I{k,t)- Gl ^ n,t 



9p,l 



(22) 



E 



p,-i 



kp\\ 



dojR J i—e,p 

dfi 



dp\\ I dpj_ 



P\\P± 
1 



<5(wr 



dfc J(M) 



9G ? ; 



\ J i 2 2 

im dp,, ^ 



ap_L o(— - W_r) 



dpn 



(23) 



7 m i 



dp,. 



with 



47rnje 2 



(17) <J(w fl 



(18) 



Performing the /c-integral in Eq. (|23j) with 
&P|| 



l + (-SL)2 



1 + (*1L)2 



-<j(fc- 



) 



(24) 



yields 



6»f 



7re 2 m, 



a [ v / i + (^) 2 



'dp. 



bill 



y =^f_ jt) : 

PI, dp,. 



Note that the summation i = e,p refers to the incoming 
beam electrons and protons. Eq. ( fl6| ) holds provided that 
the initial growth rate tp e (t — 0) due to the beam particles dGt 
is much larger than the Landau damping rate luj in the 
thermal pair plasma which is verified in Appendix B. 

We do not know precisely the wave spectrum at time 
t = 0. What we do know is that one starts with a beam 
particle spectrum that is a ^-function in momentum paral- 
lei to the ambient magnetic field and, once the instability which can be integrated over t and p f resu lting in 
stops, the particle spectrum must be a plateau in p,, . Thus 
the instability rate at late times can be effectively calcu- 
lated using the weak turbulence limit. At early times such 
is not th e case and more care has to be taken (Pohl et 
al. p002|) . However, our interest here centers on the late 
time evolution when the weak turbulence approximation 
is particularly appropriate. 



(25) 



d P\\ [G(p\\ , t) - G(p„ ,t = 0)] = irefmi 



l + (-E- f 



\i>\ 



dp 



(26) 



We then obtain 
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dt I( 



\J 1 + & 2 , dGifat') 



c7K\Y T 



E »l 







\p\ 



dp 



dp ll [G( P[l ,t)-G( Pll ,t = 0)}) (27) 



Likewise, we perform the py -integration in the wave equa- 
tion (22) using 



^ 2 m n /l + (£-)2 



wr J t—e,p 

^R.)k 

7 opj_ 



dp\\ 



dpAs^ 

7 7 TO i 



/ df l dfi 
(P-Lo P\\x — ) 



(32) 



kp\\ 



- ur) = 



-^1 + ^) 

N 2 s3/2„, m 4 c sgn(N) 



Ifcl -/V 2 - 1 



<*(P|| 



VN 2 - 1 

where N — ck/ujR is the index of refraction. We obtain 

(N 



(28) 



Eq. ( |3l| ) holds provided that the initial growth rate tp n (t = 
0) due to the beam particles is much larger than the cy- 
clotron damping rate u>i,a of the pair Alfven waves in 
the thermal background pair plasma which is verified in 
Appendix B. 

Again, we operate in the weak turbulence limit so that 
the long-time behaviour of the wave intensity can be used 
because, once again, we do not know the early time evo- 
lution of the wave spectrum. 

Transforming from the momentum variables (p±,p\i) 
to the new variables 



z.(M)-«M = o, = „ E ^T^fi 



P\\ 



E 



Pi 



1 , , ,BGi{p\\,t 
dt I(k,t)- 



dpn 



<v\\-- 



sgn(N) 



Vn 2 -i 



(29) 



rriiC 
we obtain 

4>T, 



(TTljC) 2 



(33) 



We note that precisely the left-hand side of Eq. Q27| ) taken 
at the values of p — m ^f n ^ appears in the wave equation 
(p9|), so that after inserting Eq. (|2^) we derive 



\k\ao R [^ 



uj r i=e,p 



I e {k,t)-I e (k,t = 0)=c 2 



|JV|sgn(N) 
|k|(N 2 - l) 3 / 2 



E 



L 



m i csgn(N)/v'N 2 -l 



dpi|[G(p||,t)-G(pil,t = 0)] 



(30) 



where 

_ n« 



/oo 
dE 

dfi 



dy(E 2 -l-y 2 )S(y-^+ Xl )C- dE 



•iV^)(34) 
oy 



(35) 



For Alfven speeds V e that are much less than the speed 
of light c, the index of refraction of N — c/V e >> 1 is 
large compared to unity, so that the derivative dfi/ BE is 



which relates electrostatic wave spectra and particle dis- 
tribution functions at time t to the respective quantities . . , . 

,„„„,. ,, ' r ,, negligible compared to N id h/dy), leading to 

at time t — 0. Before discussing the consequences of the 1 " ' 

quasilinear integral ([30|) for electrostatic waves we de- 
rive the corresponding quasilinear integral for transverse Wn 
waves. 



7r 2 sgn(k) 



T7 E u li( m i c f 



2.4. Time evolution of the transverse instability 



2 rd!RA_ 
R L 9wr J i—e,p 



dE 



E 



dy (E 2 -l- y 2 )6(y -- + x,)^ (36) 



Neglecting again spatial dependencies, the time- 
dependent behaviour of the intensities of the excited Transforming from the variable y to the cosine of the pitch 



transverse waves is given by (Lee 
Schlickeiser |2000|) 

din , T 
— = 2ip n I n , 

where the growth rate ip n is 

7T 



Ip |1987[ Pohl & angle p = yj\jE 2 - 1 gives 

7r 2 sgn(k) 



(31) 



,2 [9SA' 



T7 E w p,iK c ) 3 



V'n(fc) 



- .m,\v. t i E w p/' 

<9w H J i =e ,p 



dp\\ 



P 2 , 



dh 



dp± 
dfi 



el 



dEy/EP 



1 



dfi 

dp{\- p 2 )5{p- p )-^j- 



(37) 



where 



5/, 



c 1 ;u- 
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The /^-integration gives a nonvanishing value provided The two equations ( |49[ ) yield the integrals 

ImoI < 1 which is equivalent to the requirement E > E { ^ ^ ^ _ j ^ _ q^j ^ _ ^q-j 



- TV" 2 + x\ - XiN' 1 



and 



Ei l-^' 2 ~\/ 1 + x * ( 39 ) [I + (t)-I_(t)}-[I + (t = 0)-I.(t = 0)} = Z(k) (51) 

where the latter approximation holds in the limit \N\ » where 

z (*)=9iIPr E^W-V 



1. Consequently, we derive 



7r 2 sgn(k) ^ 2 ^ ^ 3 ^ ' 2|fc|c 2 r j=( 



,co ^ / #'[l-^]-^U=/*.J»(M) (52) 

^ ^V^T(l-M^)^U=, (40) n= + ,-Jo 

The influence of these excited waves on the beam particles 

In terms of the normalized phase space distribution fane- is described by the quasilinear Fokker-Planck equation for 

t lon the resonant wave-particle interaction. For Alfven waves 

/ = S(E — r)F (/i t)/[2n(m c) 3 T(r 2 — l) 1 ^ 2 ] (41) ^ e mc ^ ex °^ refraction is large compared to unity, so that 

the Lorentz force associated with the magnetic field of the 

Eq. (^) reads waves is much stronger than the Lorentz force associated 

2 with the electric field. As a consequence, on the shortest 

7rsgn(k) \ ^ ^p.i o-riju_ p-iui _..2i | t j me sca i e these waves scatter the particles in pitch an- 

zw ri au R J i= e ,p " gie /x but conserve their energy, i.e. the waves isotropise 

the beam particles. The Fokker-Planck equation for the 

wnere evolution of the phase space density is then 

^^T^^-?' ^ d -^-^\D ^1 (53) 

H denotes the step function, R t = VT/\n t \ the gyroradius. where the pitch angle Fokker-Planck coefficient is deter- 

The assumption of very large N = ck/u R -> oo, made mined by the tw0 wave intensities 
before, corresponds formally to lor — > 0, so that we can 

approximate Eq. (fl3|) as n 7T ^ i ^ Z fl } 

" mm — 2 i? 2 r 2 

=-^ = - S JEM (44) 

Mi " Mc TkV R;k 1 j r<» n 



, r , , ,. , , dkI n {k,t)S(u R - kVfi- — ) 

According to Appendix A J-00 1 

d3£A T 2c 2 fc 2 4u£ e <4 2c 2 fc 2 _ tt0 2 (1- M 2 ) f 00 Q 



fll J Wi? K - r> 2 ) 2 - u R 



( 45 ) ^ E 7rn '^ 2r2 M) /l dfcJ " (fc -^ (fcV/i + T L) (54) 



for frequencies much less than the electron gyrofrequency. . . . 

T i- /r77l\ • t-i 1 , • where we again used the limit lor — > (J. Integrating Eq. 

Inserting Eq. © m Eq. © we obtain — * * J 4 



53|) over pitch angle and time, and using Eq. ( [44[ ) we find 

*- ^ m - Rr ' 1[1 - " a f (46) f «.w.<> - ««■'.« =»)i - ^ e c- " 2 > 

1 n— H 



With the dispersion relation ( |151| ) we find for forward (+) ^ ., ( ( ,, 

(i.e. w fl = V e k) and backward (-) (i.e. u R = -V e k) Alfven dk dt I n (k,t)^6(kV^ + % = *^v\ 

waves that the wave growth rates are J-00 Jo 1 ^ 1 ' I A*l 



n=+,- ' 



and accordingly for the time evolution of the respective Eva i uat i ng Eq. © at /x = /i c gives 
wave intensities (see Eq. fey)) 

^ . +2 ^ + , ^_ _ _ 2 „_ (49) E_ I *'*.(*. on - 
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2T 2 B 2 V /" Mo i i i 
-H Mc | / dp[F i {n,t)-F i {p,t = Q)\ 



(56) 



-l 



which can be inserted into Eq. Q52| ) to yield 
V,B 2 



Z(k) 



V\k\c 2 T 



H[\k\ - R-^R 



x / dp, [Fi(fj, ,t) - Fi(n ,t = 0)] 



£ H[\k\ - RT'KiRi 



V e B 

vk 2 c 2 r 



dfi [F{n ,t)-Fi(ji ,t = 0)] 



(57) 



The system of equations ( |50| ) , (^Tj) , and ([57]) has the gen- 
eral solution at time t 



I + {t) = \jY+\(Z + I + {i))-F{Q)) 2 

+ 0.5 (Z + J+(0)-/_(Q)) 
and 

/_(*) = yy + i(^ + / + (o)-/_(o)) 2 

-0.5 (Z + /+(0)-/_(0)) 
where 

r = j+(o)/_(o) 



(58) 



(59) 



(60) 



If the initial turbulence is much weaker than the self- 
generated turbulence I(k, 0) << |Z(fc)| and has a vanish- 
ing cross-helicity l+(k,0) = I-(k, 0) = /(fc,0) we obtain 
for Eqs. (|5l,|59|) approximately 



(61) 



I ± {k,t)c\[\Z\±Z} + 



which relates the transverse wave spectra and particle dis- 
tribution functions at time t to the respective quantities 
at time t = 0. 

3. Self-excited turbulence and initial beam 
relaxation 

We first discuss which of the two, electrostatic or trans- 
verse, relaxation processes on the incoming interstellar 
electron-proton beam is faster. Hence we calculate the 
power spectra of the self-excited electrostatic and trans- 
verse turbulence from the derived quasilinear integrals for 
the beam initial condition. 

For both relaxation processes the initial boundary con- 
dition for the evolution is the same: in the beginning 
(t = 0) there is a mono-energetic beam distribution (|l|), 



i.e. in terms of the two normalised distributions and 



G Pie (pn,i = o) = a( P ||+p), 

F p , e (n,t = 0) = 6(n+i) 



(62) 



The final state of the evolution of the electrostatic insta- 
bility is reached at time t e when both growth rate and 
temporal derivative of the particle distribution function 
are zero, i.e. when dGi/dp\\ = which is referred to as 
the "plateau distribution" 



Giip^te) = ^H\p l{ + P]H[-p {] ] 



(63) 



The final state of the evolution of the transverse instability 
is reached at time t t when both growth rate and temporal 
derivative of the distribution are zero, i.e. when dFi/dfi — 
0. Consequently here 



Fi(n,t = t t ) = ~ 



(64) 



At this time the magnetohydrodynamic waves have com- 
pletely isotropised the beam distribution. 

3.1. Self-excited electrostatic turbulence 



Inserting Eqs. ( |62[ ) and (|63j) in Eq. (|3C|) allows us to inte- 
grate over pii to obtain 



/■miCSgn(N) /VN a -l 

sgn(N) / dp||[G(p||,t e )-G(p||,t = 0)] 

Jo 



H[-N] 



-.H[T 



\N\ 



Vn 2 - lVr^T Vn 2 - l 

so that the fully developed electrostatic turbulence is 

c 2 \N\H[-N] 



(65) 



I e {k,t e ) - I e (k,t = 0) 



Vr 2- ^" \k\(N 2 - 1) 



HIT- 



\N\ 1 v 

VN 2 - 1 J - 



2 2 

tot t m 2 



AirnirripC 



Vr 2 



1 



%— e,p 

m e \N\H[-N] 
m p '\k\{N 2 - If 



H[T 



\N\ 



1 



(66) 



The total enhancement in electric field fluctuation power 
due to the plateauing of the proton and electron distribu- 
tion function is obtained by integrating Eq. (|66|) 



(6E) 2 



dk[I e (k,t e )-I e {k,t = 0)] 



dN 



4irriim p c 2 m e 

~VT 2== T 7 (1 _ r - 2) - 1/2 (N 2 - l) 2 



i, r + Vr 2 - T n 



27m,- m„c 2 , m e . r „ r—z i 

; p (i + — ) rvr 2 - i - - in == 



i J 



2nTnim p c (1 
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so that the change in the electric field fluctuation energy 
density is 



AU (SE) = 



(SEl 
8tt 



-rnjm p c 2 (l + — ) 



Electrostatic waves possess equipartition between electric 
and velocity fluctuations, so that the total change in fluc- 
tuation energy density due to plateauing is 



AU = 2AU (SE) = ilXm p c 2 (l + ^) 



(69) 



which is half of the initial energy density of the beam 
particles. 

3.2. Electrostatic beam relaxation 

Using the fully developed turbulence spectrum ([36]) we 
obtain for the Fokker-Planck equation ( p5| ) in the variable 

x = \p\\\/rriiC 



at 



d 



dx Ti{x) dx 
with the time scale 



Ti(x) 



V2(r 2 - 1) 



n b 



7T(1 



-) m e m p ni p ' e xy/l + x 2 



(71) 



which is negative for all wavenumber values. According to 
Eq. ( |6l| ) this implies that predominantly backward moving 
left-handed and right-handed polarized Alfven waves are 
generated. If the initial turbulence is negligibly small the 
respective intensities of the fully developed turbulence are 



(68) I-{k,t t ) = \Z{\k\)\, l+(k,t t ) = 



(76) 



The total enhancement in magnetic field fluctuation power 
due to the isotropisation of the proton and electron distri- 
bution function is obtained by integrating Eq. (|5l]) using 
Eq. © 



/OO 
dkJ2iIn(k,tt)-In(k, 
-OO „ |_ 



0)] 



dkZ(k) = 



n=± 



m p V n^T 



47rn i rTOpV e V[l 



(77) 



so that the change in the magnetic field fluctuation energy 
density is 



(70) AU SB = 



(5Bf- 
8tt 



= \nim p V e VT[l + — 
2 rrir, 



(78) 



Alfven waves possess equipartition of wave energy densi- 
ties between magnetic and plasma velocity fluctuations, 
so that the total change in fluctuation energy density to 
pitch angle isotropisation is 



for values of < x < y/T 2 - 1. 

Evaluating this time scale at the momentum x — 
P/rriiC yields for the beam electrons 



t, 



V2 



lib 



■KUJ p ^ e to p (1 + ^) TliV 
and for the beam protons 



tp T) ( ) £g 



(72) 



(73) 



Inserting Eqs. ( |62[) and ( |64|) in Eq. 

grate over \i to obtain with 



3.3. Self-excited transverse turbulence 

d57n allows us to inte- 



(74) 



dfi [Fi(n ,t t ) 



1, 



fi(M,t = 0)] = --[l-/i c ] 



so that 



Z(k) 



2Vk 2 c 2 T 



V e B R p uj p e m 



2Vk 2 c 2 ribT m p 

(m p /m e ) \ 
R p k V 



H[\k\ - R- 1 }^^ 



(mki-R^Ki 



sgn(qj 
Rik 



Rpk 1 



H[\k\^-^R^]{l 



AU = AU&B + AU Sv = 2AU SB = n imp V e VT[l + — ](79) 

m p 

Following exactly the same calculation as in Pohl & 
Schlickeiser ( |2000| ) one can demonstrate that the increase 
( f79"l) in the energy density is balanced by a correspond- 
ing decrease in the energy density of the beam protons 
and electrons during their isotropisation. The transverse 
plasma turbulence is generated at the expense of the beam 
particles which relax to a state of lower energy density. 
The change in the energy density of the beam particles 
is given by — AU which is the fraction V e /c of the initial 
beam energy density. 

3.4. Transverse beam relaxation 

We can estimate the isotropisation length of the beam 
particles by using the fully-developed turbulence spectrum 
of backward moving transverse waves (|75|)-([76" 



J_(Mt) 

with 



i^k^HM-R- 1 ] 



V e B 2 RpLo 2 e 



2Vk 2 c 2 ribT m p 



(80) 



(81) 



(75) D 



to calculate the pitch angle Fokker-Planck coefficient (54) 
^(1-M 2 ) 



2B 2 T 2 



dk I- (k, t t )5(u>R — kV/i 



r 
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J — oo 

where 

7T Qf V e u 2 p e 



dk\k\ 



4 F 2 fipC 2 rib rrip 
The cyclotron resonance condition of beam particles of 
species i = e,p appearing in the argument of the 6- 
function in Eq. ( |S2"| ) 

lj r - Wli = Qi/T, (84) 

when combined with the pair-plasma dispersion relation 
(with V 2 <C c 2 ) for backward moving waves, 



3.4.1. Beam proton cyclotron resonance in a pair 
plasma 

For beam protons 

(83) 

7 P = m p r/m e » 1 



uj r = -V e k\ 



O 2 



r> 2 + k 2 v 2 : 



(85) 



gives the particle pitch angle cosine as a function of the 
resonant wavenumber in the form 



<>7 



ft, 



, W e + k 2 V 2 

where we introduced 

e = V e /V 



kv K r 



(86) 



(87) 

In our convention, positive [negative] k corresponds to 
right [left] handed polarization. 

It is convenient to introduce the dimensionless 
wavenumber x = k/k e and 

7* 

so that 



D — 

~ V e k 2 



(1-M 2 ) f 

J — I 



with 

g l (n,x) = ^ + (l + x 2 )- 1 / 2 + ( ll xy 
e 



(90) 



In terms of the dimensionless notation the resonance con- 
dition ( |S4|) then is equivalent to gi(/i,x) — yielding 



= -[(1 



-1/2 



(91) 



As a consequence of pitch-angle scattering the beam par- 
ticles adjust to the isotropic distribution on a length scale 
given by the scattering length 



1 -CW 



(p) 



(92) 



Eq. ( [i)2| ) is valid for scattering lengths A larger than the 
gyroradius of the particles, R p — 3 ■ lO 9 ^/!? cm and 
R e = 1.7 • 10 6 r3/i? cm for protons and electrons, respec- 
tively. If the thickness d of the outflow region is larger 
than the scattering length, indeed an isotropic distribu- 
tion of the inflowing interstellar protons and electrons with 
Lorentzfactor < T >= T(l — — T in the blast wave 
frame is efficiently generated. 

From the scattering length we derive the isotropisation 
time scale in the outflow plasma 

U = A/c (93) 



(94) 



is large and positive. According to Eq. (|9l|) for protons, 
fj. p (x) is a monotonically increasing function for x > 
going from —00 to as x goes from to 00. For x < 
0, Hp — > 00 as x — > 0—, but as —x increases, /j, p has a 
minimum at 



-7 P - 1/3 [1 



7 -2/3 r l/ 2 



X E = 

and fi p = at 

^0 = -hp - ir 1/2 - -% 

The minimum value is 

HE = -e[l-lp 2/ "? /2 ^- 



"7 



-1/3 



-e[l - ^7; 2/3 ] 



(95) 



(96) 



(97) 



Because this value is negative, there is no resonance gap 
for protons in a pair plasma even in a case of degener- 
ate cross helicity. This is in contrast with the situation 
in a cold electron-proton plasma, where one has to rely 
on electron-cyclotron waves scattering ions if waves prop- 
agating in one direction only are considered (see Vainio 
|2000 for a complete analysis of the resonance conditions 
in that case). 

The a;-value of the right-handed wave in resonance 
with the extreme value achievable by the left-handed 
waves is approximately 



da;x- 2 iJ[|a;|-e7~ 1 ]^(xg i (a;))(89) x ~ 2j~ 



1/3 



(98) 



which is always much less than unity, indicating that we 
are well in the non-dispersive range of the dispersion re- 
lation of the right-handed waves. Thus, as far as phase 
speeds are considered, we can regard all the resonant 
waves as Alfvenic. This is because dispersion modifies the 
resonance conditions enabling relativistic particles to in- 
teract with the Alfv'en waves over the whole momentum 
space. We have plotted the resonant wavenumber as a 
function of V/j, close to V [i = —V e in Fig. 1. The figure is 
plotted for protons with T = 100. 

3.4.2. Beam electron cyclotron resonance in a pair 
plasma 

For relativistic beam electrons (j e = — T e = — T < is a 
large negative number), the analysis is very similar. The 
only difference is that x has to change sign relative to the 
proton case. Now one has a monotonic resonance function, 
fj- e (x), with the left-handed waves, and the turning point 
in (J, e (x) for right-handed waves. \x e has a minimum at 

x E - r-^ii - r- 2 / 3 ]- 1 / 2 ~ r- 1 / 3 



and n e = at 

^ = (r 2 e -i)- 1/2 ^r- 1 



(99) 



(100) 



10 



Schlickeiser et al.: Conversion of relativistic pair energy 







; 

*^y *^ ' 


_ 


y 

y 

s>y y 
y^ y 


- 




y^ S 






/ / 












/ 






/ 






/ 






I 






\ 





























-2 2 4 6 

v/z + V e [V e IQ/QJ 2/3 ] 

Fig. 1. Resonant wavenumber as a function of parallel mo- 
mentum near the Alfvenic resonance gap. Here, Q = £li/T, 
v = V, the solid [dashed] curve gives k for the right- 
handed [left-handed] waves for ions and vice- versa for elec- 
trons. 



The minimum value is 



^ = - e [i-r- 2 / 3 ] 3 /^- e [i-^r- 2 /3] 



Inserting Eq. (104) into Eq. ( p2| ) gives the proton scatter- 
ing length 



X p — 



ImI 



2/3 
lp 



1 



^(l-M 2 )l(7 p 2/3 -l)^-l| 3 + 



2\ne 



iir - 1) 2 



3 V}kl 3>/2 



32 d .2/3 



C Til, ^ m p^ ^4/3 



The corresponding proton isotropisation time scale is 
A 



(105) 



+ y 

H,p — 



3\/2 1 n b iJipT 4/3 
c 8ir u> Pt e rii m e 



(106) 



Using the proton resonance approximation (|103| ) we ob- 
tain the smallest resonant wavenumbers at fi = ±1 from 
the condition g p {fi = ±1, x) = as 
c 



u p,min 



T- 



1p 



(107) 



(101) which corresponds to 



m p T v 



(108) 



Having established that no resonance gaps exist in the |fc Pjm i n | = k e \x p . m i, 
quasilinear interaction of the particles and the backward 

moving waves of relatively low wavenumbers, we will as- This indicates that the beam protons resonate with Alfven 

sume that other wave modes do not affect the dynamics of waves with wavenumbers larger than the inverse of the 

the system. In this case we can evaluate the wave spectra proton gyroradius. 
also with an alternative bi-spherical shell method which is 

described in Appendix A. „ _ . . , r , . 

3.4.4. Scattering length ot beam electrons 



3.4.3. Scattering length of beam protons 

The function gi(fi,x) appearing in the Fokker-Planck co- 
efficient (pit) for beam protons is 



g p (fi,x) = - + (1 + x )' 



-1/2 



(102) 



According to the discussion of the proton cyclotron reso- 
nance condition we approximate the function g p (fi,x) as 



9p(^, 



H + J_ for -1< a < -e 

e x~f p — f* — 



for — e < ii < 



^ + 1 + -i- for < u < 1 

e -y p x — r — 



2/3 



(103) 



allowing a straightforward approximation of the proton 
Fokker-Planck coefficient 



M 

I 2/3 



V e kl 



1)I(7p /3 



for — 1 < fi < —e 
1|~ 3 ioT-e<fi<Q (104) 
for < n < 1 



The function gi(fi,x) appearing in the Fokker-Planck co- 
efficient (B9h for beam electrons is 



+ (l + x 2 ) 



-1/2 



(T e x) 



which we approximate as 

H — for -l<u<-< 



g e (p,x) 



i 

r e x-i 

n 2/3 



r . ,_ | for -e < fi < 
H + 1 - J_ for < u < 1 



(109) 



(110) 



The electron Fokker-Planck coefficient then is approxi- 
mated as 




Vekl 



i)|(r 2/3 -i)i 



for — 1 < fi < — e 

l\~ 3 for-e<^<0 (HI) 
for < fi < 1 



Inserting Eq. (Ill) into Eq. (92) gives the electron scat- 
tering length 



A,. 



3 K 2 fc 2 = 3^2 c n bT4/3 

32 d e r 2 J 3 8?r W P,e n * 



(112) 
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and the corresponding electron isotropisation time scale 

(113) 



K 3V2 1 n b ^ i/3 



Using the electron resonance approximation (110) we ob- 
tain the smallest resonant wavenumbers at p, — ±1 from 
the condition g e (fJ, = ±1, x) = as 



±- 



which corresponds to 



— k^ 



^£ = J_ 

vT R e 



(114) 



(115) 



The beam electrons resonate with Alfven waves with 
wavenumbers larger than the inverse of the electron gy- 
roradius. Therefore, isotropisation of the beam electrons 
is possible with their own self-generated turbulence, even 
in cases where the beam protons would not generate trans- 
verse turbulence. 

3.5. Electrostatic or transverse beam realaxation: 
which is faster ? 

By comparing the respective time scales, Eqs. (f72|), (f7l|), 
( |113| ) and (106), we can answer the important question 
which beam relaxation process in pair plasmas is faster; 
plateauing in parallel momentum by the electrostatic in- 
stability or pitch-angle isotropisation by the transverse in- 
stability. 

For the beam electrons we obtain the ratio of relax- 
ation times 



t 



e,e _ " raep-7/3 << 

t t , e 3 m p 

while for the beam protons 



t 



Le ,P _ °/'"'e\l/3p-7/3 << ^ 

tt.p 3 Trip 



(116) 



(117) 



In both cases the electrostatic relaxation time is much 
smaller than the transverse relaxation time, so that the 
incoming interstellar protons and electrons first relax to 
the plateau distribution ( |63| ) in parallel momentum on 
the time sales t e ^ e and t e>v . As a consequence, the initial 
energy density of the beam particles reduces to half of its 
value; the other half has been transformed to electric field 
fluctuations of the excited electrostatic turbulence (see Eq. 

We can now study the secular evolution after time 
t e ^ e and t e ^ p , respectively, of the plateau distribution ( |63| ) 
which is still unstable with respect to the excitation of 
transverse waves. 



4. Secular transverse beam relaxation 

We again use the derived quasilinear integrals for trans- 
verse beam evolution during the secular phase. Here the 



boundary conditions are modified: at the beginning of the 
secular phase {t = t e ) the interstellar electrons and pro- 
tons obey the plateau distribution 



G i ( Plu t e ) = ^H\p ll +P}H[-p {l ] 



corresponding to 



1 



■8(fi + l)H[P-p] 



(118) 



(119) 



Again the final state of the evolution of the transverse 
instability is reached at time t s when both growth rate 
and temporal derivative of the distribution are zero, i.e. 
when dFi/dfi — 0, corresponding to 



fi(p, p,,t = t s ) 



H[P-p] 
AnPp 2 



(120) 



At this time the magnetohydrodynamic waves have com- 
pletely isotropised the plateau distribution. 

Starting from Eqs. ( |3l|) and ([lo]) we can repeat the 
earlier analysis mutatis mutandi to derive instead of Eq. 
751) that in this case 



Z(k) 



■>-\ , If 2 1 ' ■ 



\k\c 2 



'4- n? 



dpp 2 Vi 1 



P 



rriic\xi\ 



A Iff +\ ft f M 27rt^^< 

dp[h{p,p,t s ) - fi(p,fi,t e )\ = Z^TfT 



dpp 2 dp,[fi(p,fj,,t s ) - fi(p,fi,t e )} (121) 

i\Cli/k\ J-l 



Using Eqs. (|il9|)-(|l20|) we obtain 



1 - Pc 

'4tt P 2 P 



H[P-p](122) 



which after momentum integration yields for Eq. (121) 



Z(k) = - 



V e B 2 



^ 10,1 L h f. 



sgn(qi) 



2fc 2 c 2 ^ \k\Ri ' kRi 



ln(|fc|i?0] 



i 

1 sgn(qi) 



\k\Ri 



kRj 



H\k\R t )]H[\k\ - R: 



(123) 



which also results from applying the bi-spherical shell 
method of Appendix C to the flat di strib ution (119). For 
large wavenumbers \k\R p » 1 Eq. (123) is well approxi- 
mated by 



Z(k) 



V„B 2 



2k 2 c 2 VT 



(124) 



which is identical to the previous result ( |75[ ). Again Z(k) is 
negative for all wavenumber values, implying according to 
Eq. (pl|) that predominantly backward moving left-handed 
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and right-handed polarized Alfven waves are generated. If 
the transverse turbulence at time t e is negligibly small the 
respective intensities of the fully developed turbulence at 
the end of the secular evolution phase are 



I-{k,t a ) = \Z(\k\)\, l+(k,t s )=0 



(125) 



The generated wavenumber spectrum and the total change 
in transverse wave energy density is independent whether 
one starts from a monoenergetic beam (as we did in the 
initial phase of the evolution) or from a plateau distri- 
bution resulting from the faster electrostatic relaxation 
(as we did in the secular phase). This implies that we get 
exactly the same time scales for the pitch-angle isotropisa- 



tion, (106) and (113), of the plateau distribution function, 
so that at the end of the secular phase, t s = t t . e for elec- 
trons and t s = tt >p for protons, the swept-up interstellar 
protons and electrons have a flat and isotropic differential 
momentum spectrum 



N i (p,t s )=A7rp 2 f i {p,fi,t = t s ) 



H[P - p] 



P 



(126) 



in the rest frame of the plasma blob. 



5. Application to AGN jet outflows 

In order for the electrostatic and the transverse instability 
to plateau and isotropise the incoming interstellar electron 
and proton beam distributions, the respective co-moving 
time scales t e and t t have to be smaller than the crossing 
time of the beam particles which is given by the co-moving 
light crossing time of the blob 



U = d/c = 3 • 10 4 di 3 r 2 s 



(127) 



We consider these conditions separately for beam plateau- 
ing and beam isotropisation. 

5.1. Beam plateauing 



According to Eqs. ( |72| ) and (127) for beam electrons the 
requirement t e , e < U leads to 



V2 



in. 



TTU! Pi e m p n*T 3 
which is identical to 



< 3 • io 4 d:,r 2 



K, 10 ) 1/2 

n*d* 13 T 7 2 /2 



< 6.9- 10 12 



(128) 



which is well fulfilled for standard AGN jet outflow and 
environment parameters. 

Likewise with Eq. ( |73| ) for beam protons the require- 
ment t e<p < ti becomes 



K, 10 ) 1/2 



< 2.0- 10 6 



(129) 



which is also fulfilled for standard jet outflow and envi- 
ronment parameters. 



We therefore conclude that for standard AGN jet out- 
flow and environment parameters the initial beam dis- 
tributions of interstellar protons and electrons relax to 
plataeu distributions in parallel momentum, transferring 
thereby one-half of the initial energy density of the beam 
particles to electric field fluctuations of the generated elec- 
trostatic turbulence. 

5.2. Beam isotropisation and generation of radiation 



According to Eqs. (US) and (127) for plateaued electrons 
the requirement t t , e < U leads to 



Kio) 1/2 
n *d* r 7/6 

ll i u 13 L 2 



< 10 s 



(130) 



which is fulfilled for standard jet outflow and environment 
parameters. 



For plateaued protons Eqs. ( |l06j ) and (|127j) imply 



K 



b.lOJ 



a/2 



m < 4.4- 10 3 (131) 

< d 13 r 2"' 

which is also fulfilled for the adopted standard jet outflow 
and environment parameters. 

If we take the adopted AGN jet outflow (n£ 10 = 1, 
r 2 = 1, = 1) and environment (n*) parameters as face- 
value parameters, we conclude that the interstellar elec- 
trons and protons will isotropise by their self-generated 
transverse turbulence and thus be picked-up by the out- 
flow pair plasma. In Appendix D we demonstrate that 



7/6 



both transverse wave isotropisation time scales (113) and 
( |106| ) arc much shorter than the isotropisation of beam 
electrons and protons by elastic Coulomb interactions with 
pair jet plasma. 

In the pair outflow frame the external density is rii = 
Tn* and pick-up of interstellar protond and electrons oc- 
curs at a rate 



Ni(j) = wr 2 cn*y/T 2 - 15(7 - T) (132) 

The pick-up is a source of isotropic, quasi-monoergetic 
protons and electrons with Lorentz factor T in the 
jet outflow frame. By interacting with the pair back- 
ground plasma (electron bremsstrahlung, proton inverse 
bremsstrahlung), the ordered magnetic field (synchrotron 
radiation) and internal and external target photon field 
(inverse Compton scattering) the pick-up electrons and 
protons generate nonthermal photons. The detailed mod- 
elling of the radiation signature of this leptonic jet out- 
flow, including the time evolution of the energy spectra 
of radiating particles and the transformation to the ob- 
server's frame, is analogous to the earlier treatment of the 
hadronic jet outflow (see Pohl & Schlickeiser [2000 ), but 
lies beyond the scope of the present paper. We note, how- 
ever, one important difference: the pick-up protons do not 
find thermal protons in the background pair plasma to un- 
dergo inelastic p — p-collisions, so that no secondary neu- 
tral pions, secondary electrons and positrons and neutri- 
nos will be produced. These radiation products, that pro- 
vided the dominating contribution of the Dopplcr-boosted 



Schlickeiser et al.: Conversion of relativistic pair energy 



13 



GeV-TeV photon emission in the hadronic jet model and 
gave rise to correlated Doppler-boosted GeV-TeV neu- 
trino emission (Schuster et al. 2002), will be absent here. 
The pick-up protons, that carry a factor m p /m e — 1836 
more power than the pick-up electrons, can only radiate 
by the much less efficient inverse bremsstrahlung process. 
High-energy GeV-TeV emission in the leptonic jet out- 
flow model will be dominated by inverse Compton scat- 
tering of pick-up electrons off internal and external tar- 
get photon fields. It will be most interesting to study in 
future work quantitatively the leptonic jet radiatio n sig - 
natures modifying the formalism of Bottcher et al. ( [2001 ) 
developed for optically-thick collision-dominated relativis- 
tic pair gamma-ray burst fireballs. 

5.3. Build-up of transverse magnetohydrodynamic 
turbulence 



The co-moving time scales 



0.19 



b.lOJ 



a/2 



1/6 



and 



^) 4/ % e = 4.3 • 10 : 



(n* ) 1/2 

, V L b,w) 

„*r ,1 / e 



(133) 



(134) 



given in Eqs. (113) and (106), respectively, also provide es- 
timates for the time to build up the fully developed power 
spectrum of transverse hydromagnctic turbulence from 
the beam electrons and the beam protons, respectively. 
Hydromagnetic turbulence upstream and downstream is 
crucial for diffusive shock acceleration to operate. We re- 
call that beam electrons alone only produce transverse 
turbulence at wavenumbers larger than the inverse of the 
electron gyroradius (|&| > R^ 1 ) with an (integrated) en- 
ergy density smaller by a factor (m e /m p ) — 1/1836 than 
the turbulence produced by the beam protons. As our ear- 



lier estimate (13C) indicated, this small fraction of electron 
generated transverse turbulence is built up after time t ti e 
for standard AGN outflow and environment parameters. 
For the protons to provide the much larger fraction of 
transverse turbulence, the longer time t tjP is needed, and 
it will occur at later ti mes of the AGN jet evolution. The 
co-moving time scales (133)-(134) correspond to Doppler- 
shortened time scales t* t =t t / Dl , where the Doppler fac- 
tor D L = T~ l [l - (Vcos0*/c)] -1 , depends on the ob- 
server's viewing angle 9* . 

Eqs. ( |l33| ) and ( |l34[ ) indicate that it takes a non- 
negligible fraction of the early AGN jet phase, to build 
up the magnetohydrodynamic turbulence for nonthermal 
particle acceleration at either internal or external shocks 
associated with the pair jet outflow. 

6. Summary and conclusions 

We have investigated the microphysical details of the en- 
ergy conversion in relativistic pair outflows interacting 



with the surrounding interstellar medium consisting of 
cold protons and electrons. We have represented the rel- 
ativistic pair blast wave as a one-dimensional channeled 
outflow directed parallel to a uniform background mag- 
netic field. Viewed from the coordinate system comov- 
ing with the pair outflow, the interstellar protons and 
electrons represent a proton-electron beam propagating 
with relativistic speed antiparallel to the uniform mag- 
netic field direction. We demonstrate that the beam ex- 
cites both electrostatic and low-frequency magnetohydro- 
dynamic Alfven-type waves via a two-stream instability 
in the pair background plasma, and we calculate the time 
evolution of the distribution functions of the beam par- 
ticles and the generated plasma wave turbulence power 
spectra. 

For standard AGN jet outflow and environment pa- 
rameters we show that the initial beam distributions of 
interstellar protons and electrons quickly relax to plateau- 
distributions in parallel momentum, transferring thereby 
one-half of the initial energy density of the beam particles 
to electric field fluctuations of the generated electrostatic 
turbulence. 

On considerably longer time scales, the plateaued in- 
terstellar electrons and protons will isotropise by their 
self-generated transverse turbulence and thus be picked- 
up in the outflow pair plasma. These longer time scale are 
also characteristic for establishing fully developed power 
spectra of transverse hydromagnetic turbulence from the 
plateaued electrons and protons. This hydromagnetic tur- 
bulence upstream and downstream is crucial for diffusive 
shock acceleration to operate at external or internal shocks 
associated with pair outflows. It takes a finite time period 
of the early AGN outflow phase to build up the magneto- 
hydrodynamic turbulence for nonthermal particle accel- 
eration at either internal or external shocks associated 
with pair jet outflows. During this initial phase - prior 
to the generation of hydromagnetic turbulence - the ra- 
diative and energy-exchange processes in the jet outflow 
are not dominated by acceleration processes of nonthermal 
charged particles at collisionless shock waves. 
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7. Appendix A: Dispersion relations of parallel 
plasma waves in a thermal magnetised pair 
plasma 

The properties of parallel propagating plasma waves 
in magnetised thermal plasmas have been discussed by 
Fichtn er & Schlickeiser ( 1995 ) and Schlickeiser et al. 
( |l997| ). For subluminal waves with index of refraction 
N = ck/ujR > 1 the dispersion relations of longitudinal 
(Al = 0) and transverse (At = 0) can be split in their 
real and imaginary parts. 

7.1. Longitudinal waves 

According to Eq. (22) of Fichtner & Schlickeiser ( |1995[ ) we 
obtain for longitudinal waves in the limit of weak damping 

UJ = UJr + IUJJ, UJj « UJr, 



ftA L (uj R ,k) = 1 



a 2 {l + 2a 2 



dy(l 



L0 2 R K 2 (Ha) Jo 



and 



3A l (cjr,/c) = ~sgn(w R ) 



V) 



K 3 (» ay /l + cT2-y* 



(1 



y 2)3/2 

(135) 



la ex P(~faVl + cr 2 ) 
c4N 3 Ai a K 2 (/i a ) 



[2 + 2^ a Vl + a 2 + ^ 2 a (l + a 2 )] =0 



(136) 



where a — (N — l)" 1 / 2 , /i a = 1/0 O = m a c 2 /(k B T a ), and 
uj Pj a de notes the plasma frequency of species a. Solutions 
of Eq. ( |135] ) define the frequency-wavenumber relations 
= oJn(k) of all possible modes in the system whereas 
the growth or damping rate of each mode is calculated 
from Eqs. (135) and (136) by the standard weak-damping 
expansion 



^sK L (uj Rl k) 



(137) 



The integral equation ((135) has been further reduced 
by Schlickeiser & Mause ( |l997 ). Introducing normalised 
wavenumbers k — ck/u) Pte and frequencies / = ojR/oj p . e 
these authors obtained for non-relativistic temperatures 
Ha » 1 (see Eqs. (34)-(38) of Schlickeiser & Mause 
(Il997|) ) the approximation 



RA l (/,k)~1- 



3^ 2 



2/ia Haf 2 



(138) 



with £ a = cJp )0 /cJ P)e . For a single temperature plasma this 
approximation holds for the wavenumber range 



[y>+^)] i/2 <M<2/4 /2 [i ■' 



a 



16/i 



For a single-temperature (/x a = l/Q) non-relativistic 
electron-positron plasma to lowest order in the normalised 
temperature 9 << 1, Eq. (138) reduces to 



KA L ( Wit ,fc)~l 



yielding electrostatic waves 



(140) 



(141) 



in the wavenumber range |A;| < 2w Pie /(c9 1 / 2 ) as only lon- 
gitudinal plasma mode. Moreover, we obtain 



6®lk L {u>R,k) 
duj R 



2 



(142) 



While superluminal electrostatic waves are not damped, 
in the subluminal wavenumber range V2uj p . e /c < \k\ < 
2t£ Ei e/(c0 1 / 2 ) we derive for the associated damping rate 
(II 



-7rsgn(w R ) 



exp(-Vl + CT 2 /e)[2 + 
= -7rsgn(w R ) 



2yT 



c 3 k 3 K 2 (l/6) 



e 



e 2 



U>1 



exp 



ck[c 2 k 2 - 2w 2 je ]6K 2 (l/e) 
c\k\ 



26 2 (c 2 fc 2 -2u; 2 e ) 



c 2 k 2 



(143) 



7.2. Transverse waves 



The dispersion relation of parallel propagating transverse 
waves in a thermal plasma with non-relativistic temper- 
atures has been analysed in Sect. 3 of Schlickeiser et al. 



( |1997D . Its real part 
5?A T 



1 



N 2 



E 



kv t h,aVR 



lev 







th,< 



(144) 



involves the plasma dispersion function Z(x) (Fried & 
Conte 1962 ), the non-relativistic gyrofrequencies Sl a — 
q a B /(m a c) and the non-relativistic thermal velocities 
v th,a = 2c 2 //x a = 2k B Ta/m a . In Eq. ( p.44|) we adopt the 
convention that positive values of ljr > denote left- 
handed circularly polarized waves, and that negative val- 
ues of ojr < denote right-handed circularly polarized 
waves. An expansion of Z(x) for large arguments and the 
weak-damping approximation \ujj/ujr\ << 1 allows the 
derivation of a simplified, asymptotic dispersion relation 
for transverse oscillations propagating parallel to the ex- 
ternal magnetic field 



(139) 3?A T ^ 1 



Ur(ur - Q a ) 



= 



(145) 
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According to Eq. (54) of Schlickeiser et al. (1997) we ob- 
tain for the weak damping rate in the same limit 



tt 1 / 2 

WR 



-sgn(w R ) 



a?RA T (wR,k) 



exp( 



ai2\ 



(146) 



For a single-temperature (/x a = 1/6) non-relativistic 
electron-positron plasma Eq. (145) reduces to 



2 1.2 



5RA T ^ 1 



c 2 fc 



-02 







(147) 



For subluminal phase speeds (\ujR/k\ << c) we may ne- 
glect the first term on the right hand side of Eq. (147), so 
that the transverse dispersion relation yields 



fi 2 + K 2 fc 2 



(148) 



where 



V, = 



fi,. 



= B [87rnf,m e ] 



-1/2 



6.6 ■ 10^( J B/Gauss)(n b /cm-'7 
9.2 • 10 7 (e B ^ 1 e_ 4 ) 1/2 cms" 1 



- 1 / 2 cm s- 1 



(149) 



is the Alfven speed in the pair plasma. For wavenumbers 
much less than |fc| << k e , where 



k e — fi e /Ve 



2.7 • W^n 1 ^ cm" 



(150) 



is the inverse pair skin length, Eq. (148) reduces to 



= V 2 k 



21.2 



'R-v*» (151) 

whereas in the opposite case \k\ » k e Eq. (|14§) ap- 
proaches the pair-c yclo tron wave limit u> R ~ VL 2 . 
Moreover, Eq. ( |147|) yields 



dUA T (to R ,k) 4wl e w R 2c 2 k 2 



dujR 



+ 



(152) 



J R 



so that the damping rate ( 146 ) becomes 

/; \ 1 r~n~ Up sgn(w R ) 
W/(fc) = ~2 V 26^7^ 



exp( 



(tJfl - fi e ) 

26fc 2 c 2 



exp( 



(lor + fi e ) 2 
26fc 2 c 2 



(153) 



Especially, for small frequencies (\u>r\ << fi e ) we obtain 
for the pair Alfven waves (|151[) the damping rate 



rV 2 



fi! 



fi e 



2 «t M |A|(l 



k 2 

k 2 



Vth,e 



(154) 



where t> t /, e = v2Bc = -\/2fcsT e /m e is the non-relativistic 
thermal pair velocity. 

To check the weak-damping approximation we calcu- 
late 



TT 1 / 2 , P 



fc 2 J 



(155) 



where we introduce the pair plasma beta (3 P = V 2 /v 2 h = 
0.05es._i. One notices that for small wavenumber values 
|fc| << [3p^ 2 k e = 0.22e^/ 2 _ 1 fc e the weak-damping approxi- 
mation is well satisfied. 



8. Appendix B: Comparison of the initial growth 
rates with the thermal Landau and cyclotron 
damping rates 

We noted that plasma waves, due to the unstable inter- 
stellar beam particles entering the jet outflow pair plasma 
will only grow if the initial growth rate is larger than the 
thermal Landau (143) and cyclotron (154) damping rates. 
Here we investigate under which conditions this is the case. 
We consider longitudinal and transverse waves in turn. 

8.1. Electrostatic waves 

If we include the contributions from the beam particles 
the initial longitudinal dispersion relation reads 



A e (t = 0) = 3f?A L + i^Al 



2tt 

10 



p. i 



dpu 



dpA 



P±P\\ 



j(u> — kvu) 



i 

dfi(p.L,P\\,t = 0) 
dpii 



(156) 



where !RAl and 3Al are given by Eqs. (14C) and (136), 
respectively. Neglecting 3A^ for this calculation we obtain 
in terms of the normalised phase space distribution ( pl| ) 



A e (t = 0) ~ 9?A L + 



dp\\ p\\ 




kp\\ 



dGj(p h t = 0) 
dp\\ 



(157) 



Inserting the initial distribution function (|6^), i.e. 
G(p\\,t = 0) = S(jpu + P), and integrating partially with 
respect to p» we derive with Eq. (140) 



Ae(f = 0) = 1 - 



CO 2 



p,i 



r 3 (w + kvy- 



= o 



With the abbreviations 



= y/2u 



/2 



and the substitution 
/ = to + kV 



(158) 



(159) 



(160) 
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Eq. (15S) becomes 
\kV-f\ 



(161) 



2e 2e 2 n tt 1 / 2 . . r 
— + T2-] - -^e37i^e r ex P[ — 



e 



(171) 



A posteriori it can be shown that in the resonant wave 
regime 

w 2 « T 3 / 2 (162) 



For the ratio of the maximum growth rate to this damping 
rate we obtain 



r - 1, 



\ui{k )\ 



= v/3M^) 1/3 e 3 / 2 r- 7 /3 exp( „ — ) = 



'6,10 



e 



so that the sqare root in Eq. (161) can be expanded: 

,,2 



«) 1/3 Kio) 1/3 0-4 r 2 7/3 lO 4 - 23 - 1O5r2e -- 18 » 1 (172) 

which is a huge factor for standard AGN jet outflow 
and environment parameters due to the exponential de- 
pendence on exp[(r — l)/0]. So we can safely neglect 



kV ~ / - w e - 



2r 3 / 2 



(163) 



This equation has already been examined by Buneman ,, T , , . c , / , , , , , , . 

,ii „ , . . the Landau damping of the beam generated electrostatic 

( |1959| ). Substituting 



1/3 2/3 
p W e W^ 

f = x, W 



-W e 



1/3 2/3 
We Wi y 



Eq. ( |163| ) reduces to 
x 3 - yx 2 - 1/2 = 



waves. Only for values of < 4 • 1O~ 5 0_4, i.e. very 
hot pair plasmas and slow Lorentz factors of the outflow, 

(164) the damping of electrostatic waves outnumbers the initial 
wave growth. 

(165) 8.2. Transverse waves 



This cubic equation has one real and two complex conju- Including the contributions from the beam particles to 
gate solutions x\^ provided y > —3/2, corresponding to the initial transverse dispersion relation we obtain after 
wavenumbers 



k > 



o 1/3 2/3 
Sw e «V 

[ 2TV 



W e 

~v 



The complex conjugated solutions are 



y 1 

Xl < 2 =3~ 2 



.1 



.1 



ft + __ + A 1 ^)!/ 3 + (- + — - A 1 / 2 ) 1 / 3 
V 4 27 ' V 4 27 ' 



±i 



\/3rA 



2 

where 



^4 + 27 



A X /2)l/3 



(i + l-A 1 / 2 )^] ( i 6 7) 



1 y 

A = h — 

16 54 

The maximum growth rate Sav, 
y = 0, corresponding to 

fc = -w e /V = -2 l ' 2 u> p , e /V 

Consequently 



(168) 

= 3 1 / 2 2- 4 / 3 occurs for 
(169) 



1/3 2/3 
We W,' , 



3I/2 



r 



1/3 



2r 



- [ ^ ( i + ^ ) ]i/3 = ^! Wpie[ ^L ] i/3 

1 lib m p 2 Lni 



(170) 



analogous manipulations as in the previous subsection 

- 2 fc 2 2^ 2 e UjleM^ + kV) 



(166) A t (t = 0) = l 



r me 



1 



kV 



r 



ui 2 Tni> 

-■ (173) 



where the two beam contributions stem from the interstel- 
lar protons and electrons, respectively. The beam protons 
will excite transverse waves near their resonance 



kV - 



TTle^le 



= 



(174) 



whereas the beam electrons excite transverse waves near 

fie 

r 



kV • 







(175) 



Here we consider only the waves generated by the inter- 
stellar protons, and leave the case of electron generated 
waves to the interested reader. We also assume small wave 
frequencies, i.e. uo 2 « f2 2 so that Eq. (173) reduces to 



A t (t = 0) ~ 1 



.21.2 



C 2 k 



C 

V 2 



ul,e n * 



kV 



a; 2 r?ib 1 



kV 



(176) 



Now we compare this maximum growth rate with the This approximation holds for wavenumbers fc near the pro- 



Landau damping rate (143) calculated at the wavenum- 
ber kg which is 



ton resonace condition (174). Using V 2 « c 2 we then 
obtain 



tt (r 2 

23/2 W P> e 



1)3 /2 



-rye 



QK 2 {l/Q) 



[1+ 



= lj 2 - V 2 k 2 - 



kV 



p,e 



kV 



(177) 
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where we introduce the parameter 

Q = vim 

c 2 Tni, 



Eq. (177) leads to the cubic equation 

oj 3 + buj 2 + gco + h = 

with 



(178) 



(179) 



b = kV - O, g = 
h = kVg + nV?k 2 , = 



(V 2 k 2 + au; 2 J, 



Provided 

D = Q 2 + S 3 > 0, 

where 

Q = (b/3) 3 - (bg/6) + (h/2) = 
l {kV _ 0)3 + 2 _ kVg _ Qg + \w 2 k 2 



(180) 



(181) 



and 

S=(g/3)-(b/3) 2 



[(fcV-f2) 2 + 3.g], 



(182) 



(183) 



the cubic equation ( |179| ) has two complex conjugate solu- 
tion 

u = \(Sl-kV)- h(-Q + V^) 1/3 + (-Q - ^D) 1/3 \ 



(184) 



The maximum growth rate is a measure for the impor- 
tance of the instability. For its estimate we use the ansatz 



of Achatz et al. (1990) that it is largest when the reso- 



nance condition (174) is fulfilled yielding with lor = ±V e k 
for relativistic outflow velocities (V >> V e ) 



k 



n 
v 



V±V e 

For this wavenumber we obtain 



(185) 



S(k ) = ~au£ e [l + 3r], Q(k ) = ^au£ e [l + (186) 

with the frequency ratio 

V 2 fl 2 2m 2 , V 2 nt 



aul e V 2 m 2 p c 2 T(T 2 - 1) n* 



5. 6 -io- 8 ^e S) _ 1 e_ 4 r 2 - 3 



(187) 



which for standard jet outflow and environment parame- 
ters is very small. We then obtain 



[(1 



Q(fco) 



a/3 



(1- 



Q( k o) U/3j 



D^ 2 (k )- 



3ft 2 



^/ 2 D 1 / e (k )^(a/3) 1 ^[l + ^} 1 / 6 



because 

Q 2 (k ) 
D(k ) ' 



(Ve/vy 



< 1 



(188) 



(189) 



is always smaller than unity. 
Collecting terms we obtain 



oJi(k ) 



V3 c n* 



1-4 f2 e ,.n*IV 1/2 r-, , n „,„ n U0 r 2- 



io- 4 ^( 



•[1 + 0.018- 



(190) 



6,10 



which has to be compared with the cyclotron damping 
rate (154) calculated at the wavenumber kg: 



TT 1 / 2 ,2 

i^i,A(k ) = — tt e die 1 



(191) 



where d\ = 1.3-10^26— 4 -1 / 2 . Because of the small factor 
di exp(— d\) = exp(—1.7-10 14 ) the cyclotron damping rtae 
is negligibly small compared to the initial growth rate. 

9. Appendix C: Calculation of the wave spectrum 
with the bi-spherical shell method 

Here the self-generated wave spectrum of backward mov- 
ing transverse waves can be evaluated with the assump- 
tion, that the free energy (as measured in the plasma rest 
frame) of the initial beam particle distribution is given 
to these waves in the isotropization process. This method 
has been used before for non-relativistic particles when 
evaluating the wave spectrum generated by freshly ionized 
pick-up ions (e.g., Huddleston & Johnstone 1992; Isenberg 
& Lee 1996). Generalizing this method to relativistic par- 
ticles and waves with non-relativistic phase speeds goes 
straight-forwardly by replacing particle speeds with their 
momenta. 

The final particle distribution is uniform in the so- 



called (Huddleston & Johnstone 1992) bi-spherical shell, 
which in our case of almost Alfvenic waves and particles 
injected at p± = is accurately given by the spherical 
shell 



2,2 2 
P± +P\\ =P , 



(192) 



with the momentum measured in the Alfven wave frame, 
where its magnitude, p, is taken as constant. (The final 
distribution function is, thus, isotropic in the Alfven-wave 
frame.) We may write the self-generated wave spectrum 
in form (Isenberg & Lee 1996) 



dE 

I a (k) = 4n^-N a {p\\(k)}, 



(193) 
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Fig. 2. Spectrum of self-generated waves for after a rela- 
tivistic pick up of electrons and protons with a T = 100. 
The left [right] panel gives the spectrum of the right- 
handed [left-handed] waves. Electron-generated (proton- 
generated) waves are plotted with a dashed curve (solid 
curve) . 

where dE/dk = \dE / dp^\\dp^ / dk\ is the energy lost by 
a particle (per unit resonant wavenumber) as it scat- 
ters through the parallel momentum range resonant with 
wavenumber range (fc, k + dk), and 



N a {p\\(k)} = n a 



P-P\\(k) 
2p 



(194) 



is the number density of particles scattered accross the 
parallel momentum p\\(k), resonant with k, during the 
isotropization. We have \dE/dp\\ \ = V e (with E still mea- 
sured in the plasma frame, of course), so 



Ia(k) = 2-Kn a V e 



dp 1 1 



dk 



P~P\\ 
P 



where one uses Eq.(|8q) (and py 
Alfven-wave frame to get 



P\\(k) 



k 



2k 2 



(195) 

PHi) evaluated in the 
(196) 



for the fc-rangc of interest to our analysis. Thus, the spec- 
trum reads 



I a (k) = 2Tm a pV e 



, V e k 



1 



R L , a k 2 



RL, a k' 



(197) 



where Rl.u — lv/\Q, a \ is the particle's Larmor radius and 
e Q = sign(il Q ). The small correction due to the dispersive 
effects in dp\\ j dk has to be retained, but the corresponding 
term in p — py can be neglected, since it is small in the 
region p» p. 

The above spectrum holds only for a single particle 
species being picked up by the waves. When several par- 
ticle species are picked up simultanously, then it may not 



be possible to avoid considering multiple resonances. For 
the least massive species the reasoning behind the neglect- 
ing of them is still valid, but for the more massive ones, 
part of the free energy of the lightest species will go to the 
"denied" bands of the massive particles, which may inval- 
idate the assumption of single resonant wavenumber. In 
the specific case of electrons and protons, however, the re- 
gion where multiple resonance are possible for protons, the 
electron-generated spectrum has a steep slope (oc k~ 2 ), 
which as a first approximation still enables us to neglect 
multiple resonances in case of protons, as well. Thus, the 
total wave spectrum comes from a contribution from pro- 
tons and electrons as 



I(k) =2irY, H (\ k \RL, a -l)n aPa V e 



1 



V e k w 



R L , a k 2 



(198) 



Taking n a = m, v a — V, p a = Tm a V with T = 
- {V/c) 2 , and Rh.ct = TV/\Cl a \, gives 



Ak e rib 



k 2 e Tm a k \ ( n a \ 



k 2 " m e k e 



(199) 



showing that the electronic contribution at |O e |/rV|fc| <SC 
(me/TmpY^ke -C fce/r 1 / 3 is of the same order as the 
ionic one. The spectrum and its constituents are plotted 
in Fig. |. 

10. Appendix D: Isotropisation by Coulomb 
interactions 

Apart from scattering with transverse plasma waves addi- 
tional pitch angle scattering of the platcaued interstellar 
electrons and protons is provided by the elastic Coulomb 
interactions with the pair electrons and positrons in the jet 
outflow plasma. Here we demonstrate that the Coulomb 
isotropisation time scales of the beam electrons and pro- 
tons are much longer than the plasma wave isotropisation 

time scales (106) and (113). 

According to Kulsrud et al. ( 1993 ) the linearized 



Fokker-Planck scattering operator for Coulomb interac- 
tions is 



(200) 



where the electron scattering rate due to M0ller and 
Bhabha scattering (Bottchcr et al. 2001 ) is 



= 2.6 • 10~ 3 



'b.lO 1 - 2 



r- 1 Hz 



(201) 



while the proton scattering rate (Mannheim & Schlickeiser 
1994] ) is 

v ep = 3.3 ■ 10~ 8 KioiV 1 Hz ( 202 ) 
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The solution of the associated Fokker-Planck equation 
during the secular phase 



(f ) - (§) 



at > c 



with the initial condition (11£), 
1 



fi(P,fJ-,t e 



^ + 1)H[P- V ] 



(203) 



(204) 



2irPp 2 

is easily obtained in terms of Legrendre polynomials as 
H[P-p] 



fi(P,fJ-,t > t e ) 



AnPp 2 



x ]T(-l) n (2n + l)P n (p) exp[-n(n + \)v el (t - t e )] (205) 

n=Q 

For times larger than 

tc,i = {Ive^- 1 (206) 
Eq. ( |205| ) approaches the isotropic distribution ( |120| ) 
H[P-p] 



hip) 



AirPp 2 



(207) 



here without generating transverse waves. 

For electrons the Coulomb isotropisation time scale is 



t c ,e = 1-9 -10 2 T 2 (n^o)" 1 
and for protons 
t CiP = 1.5 • 10 7 T 2 (n 



(208) 



(209) 



We note that for both partic les t hese ti me s cales are much 
longer than the time scales (|ll3| ) and ( |l06|) for isotropisa- 
tion by transverse plasma waves calculated for standard 
AGN jet outflow and environment parameter 



and 



3.0-10" 



(<io) 



1/2 



^ r 2 



S. 



( n t io) 1/2 

k,p = {m P /m e ) A/:i tt,e = 6.7 ' , ,„ s 



n*T 



1/6 



(210) 



(211) 



Dermer, C. D. & Schlickeiser, R., 1993, ApJ 416, 458 
Dermer, C. D., Sturner, S. J. & Schlickeiser, R., 1997, ApJS 
109, 103 

Fichtner, H. & Schlickeiser, R., 1995, Phys. Plasmas 2, 1063 
Fried, B. D. & Conte, S. D., 1961, The Plasma Dispersion 

Function, Academic Press, San Diego 
Hartman, R. C, et al., 1999, ApJS 123, 79 
Huddleston, D. E. & Johnstone, A. D. 1992, J. Geophys. Res. 

97, 12 217 

Isenberg, P. A. & Lee, M. A. 1996, J. Geophys. Res. 101, 11 055 
Kulsrud, R. M., Sun, Y.-G, Winsor, N. K. & Fallon, H. A., 

1973, Phys. Rev. Lett. 31, 690 
Lee, M. A. & Ip, W.-H., 1997, J. Geophys. Res. 92, 11041 
Lerche, I., Pohl, M. & Schlickeiser, R., 2000, J. Plasma Phys. 

64, 459 

Mannheim, K., 1993, A& A 269, 67 

Mannheim, K. & Biermann, P. L., 1992, A& A 253, L21 
Mannheim, K. & Schlickeiser, R., 1994, A& A 286, 983 
Maraschi, L., Ghisellini, G. & Celotti, A., 1992, ApJ 397, L5 
Pohl, M., 2001, Rapporteur pape r at the 27th Intern. Cosmic 

Ray Conf. ( |astro-ph/0111552| ) 
Pohl, M., Lerche, I. & Schlickeiser, R., 2002, A& A 383, 309 
Pohl, M. & Schlickeiser, R., 2000, A& A 354, 395 
Schlickeiser, R., 2002, Cosmic Ray Astrophysics, Springer- 

Verlag, Heidelberg 
Schlickeiser, R., Fichtner, H. & Kneller, M., 1997, J. Geophys. 

Res. 102, 4725 

Schlickeiser, R. & Mause, H., 1995, Phys, Plasmas 2, 4025 102, 
4725 

Schuster, G, Pohl, M. & Schlickeiser, R., 2002, A& A 382, 829 
Sikora, M., Begelman, M. C. & Rees, M. J., 1994, ApJ 421, 
153 

Sikora, M. & Madejski, G, 2001 in: Current Hig h Energy 
Emission around Black Holes (astro-ph/0112231 ) 

Tavecchio, F., Maraschi, L. & Ghisellini, G, 1998, ApJ 509, 
608 

Vainio, R. 2000, ApJS 131, 519 

Vainio, R. & Schlickeiser, R., 1999, A& A 343, 303 

Vainio, R. & Schlickeiser, R., 2001, A& A 378, 309 



References 

Achatz, U., Lesch, H. & Schlickeiser, R., 1990, A& A 233, 391 
Arbeiter, G, Pohl, M. & Schlickeiser, R., 2002, A& A 386, 415 
Atoyan, A. & Dermer, C. D., 2001, Phys. Rev. Lett. 87, 221102 
Bednarek, W. & Protheroe, R. J., 1999, MNRAS 302, 373 
Bloom, S. D. & Marscher, A. P., 1996, ApJ 461, 657 ApJ 563, 
71 

Bottcher, M., 2001, in: Gamma-Ray Astrophysics through 
Multiwavelengt h Experiments 20 01, Bull. Astron. Soc, 



India, in press (astro-ph/0105554) 



Bottcher, M., Mause, H. & Schlickeiser, R., 1997, A& A 324, 
395 

Bottcher, M., Schlickeiser, R. & Marra, A., 2001, ApJ 563, 71 
Buneman, O., 1959, Phys. Rev. 115, 503 



